arXiv:1504.01624v3 [hep-ph] 10Jul2015 


Eur. Phys. J. C manuscript No. 

(will be inserted by the editor) 


Solving the inhomogeneous Bethe-Salpeter Equation in 
Minkowski space: the zero-energy limit 


Tobias Prederico Giovanni Salme El2l^ Michele Viviani ISEl 

^Dep. de Ffsica, Instituto Tecnologico de Aeronautica, DCTA, 12.228-900 Sao Jose dos Campos, Sao Paulo, Brazil 
^Istituto Nazionale di Fisica Nucleare, Sezione di Roma, P.le A. Moro 2, 1-00185 Roma, Italy 
^Istituto Nazionale di Fisica Nucleare, Sezione di Pisa, Largo Pontecorvo 3, 56100, Pisa, Italy 


Received: date / Accepted: date 


Abstract For the first time, the inhomogeneous Bethe- 
Salpeter Equation for an interacting system, composed 
by two massive scalars exchanging a massive scalar, is 
numerically investigated in ladder approximation, di¬ 
rectly in Minkowski space, by using an approach based 
on the Nakanishi integral representation. In this pa¬ 
per, the limiting case of zero-energy states is considered, 
extending the approach successfully applied to bound 
states. The numerical values of scattering lengths, are 
calculated for several values of the Yukawa coupling 
constant, by using two different integral equations that 
stem within the Nakanishi framework. Those low-energy 
observables are compared with (i) the analogous quan¬ 
tities recently obtained in literature, within a totally 
different framework and (ii) the non relativistic evalu¬ 
ations, for illustrating the relevance of a non perturba¬ 
tive, genuine field theoretical treatment in Minkowski 
space, even in the low-energy regime. Moreover, dy¬ 
namical functions, like the Nakanishi weight functions 
and the distorted part of the zero-energy Light-front 
wave functions are also presented. Interestingly, a highly 
non trivial issue related to the abrupt change in the 
width of the support of the Nakanishi weight function, 
when the zero-energy limit is approached, is elucidated, 
ensuring a sound basis to the forthcoming evaluation 
of phase-shifts. 

Keywords Bethe-Salpeter equation • Minkowski 
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1 Introduction 

Within a field theoretical framework, it is a highly non 
trivial challenge to develop non perturbative tools in 
Minkowski space, but it is quite desirable to devote ef¬ 
forts in that direction, in order to gain insights that 
could turn out useful in particle physics. In the last few 
years, solving the homogeneous Bethe-Salpeter equa¬ 
tion (BSE) |T], directly in Minkowski space, has made a 
substantial step forward pinTlHlIHlIHirTH^lUi nTillTT] due to 
approaches based on the so-called Nakanishi perturbation- 
theory integral representation (PTIR) of the n-leg tran¬ 
sition amplitudes [12) . 

The Nakanishi PTIR for the three-leg amplitude is 
emerging as a very effective tool for studying the bound 
state problem [2l[8ll^[5ll6l[7l[8l [TT)ll 11) . within a rigorous 
field-theory framework. Though the Nakanishi PTIR of 
the three-leg amplitude, or vertex function, had been 
devised within the perturbative framework of the Feyn¬ 
man diagrams (as it happens for any n-leg amplitude 
PTIR), it has been shown to work extremely well as 
the initial Ansatz for obtaining aetual solutions of the 
homogeneous BSE. It must be recalled that BSE, be¬ 
ing an integral equation, belongs to a non perturbative 
realm, and therefore the Nakanishi integral representa¬ 
tion of the three-leg amplitude can be only an Ansatz, 
when exploited in this context. 

The main features of the Nakanishi integral repre¬ 
sentation of any n-leg amplitude are basically related 
to the formal infinite sum of the parametric Feynman 
diagrams, that contribute to the amplitude under con¬ 
sideration. In particular, the n-leg amplitude PTIR has 
a well-defined structure, given by the folding of (i) a 
denominator, containing all the allowed independent 
invariants and governing the analytic behavior of the 
amplitude itself, and (ii) a weight function, that is a 
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real function depending upon real variables (one is a 
non compact variable, while the others are compact). 
It should be emphasized that, at this stage, the Nakan- 
ishi weight function has only a formal expression HU. If 
there were an equation for explicitly determining such a 
weight function, then one could quantitatively evaluate 
the actual n-leg amplitude, under consideration. The 
homogeneous BSE, that obviously does not belong to 
the original framework of PTIR, has inspired a differ¬ 
ent usage of the formal expression of a particular n-leg 
amplitude, namely the three-leg one, or vertex func¬ 
tion. Indeed, if one assumes that the Nakanishi integral 
representation of the three-leg amplitude be formally 
valid also for the BS amplitude (still a three-leg am¬ 
plitude, but for a bound state), then the weight func¬ 
tion could be considered as an unknown function to 
be determined. It has to be pointed out that, a pri¬ 
ori^ there is no guarantee that such an approach for 
solving BSE be successful, given the caveat above men¬ 
tioned. Fortunately, it works, as shown in Refs. Ham 
[aiHieiiniiii], where the above strategy was applied, 
but with some differences, for solving the homogeneous 
BSE directly in Minkowski space. More precisely, by 
using the PTIR Ansatz for the BS amplitude one can 
derive, in a formally exact way, an equation for the 
Nakanishi weight function, starting from the homoge¬ 
neous BSE, and look for solutions. If the new equa¬ 
tion for the weight function has solution, then one can 
claim that BS amplitudes, actual solutions of the ho¬ 
mogeneous BSE in Minkowski space, can be (i) for¬ 
mally written like the PTIR three-leg amplitude, and 
(ii) numerically determined. In order to achieve a for¬ 
mally exact integral equation for the weight function 
from BSE, it is very useful and effective adopting a 
Light-front (LF) framework. This has been done both 
in the covariant version of the LF framework H] ^-nd 
in the non-explicitly covariant one [9j. In particular, 
the bound states of a massive two-scalar system inter¬ 
acting through the exchange of a massive scalar have 
been studied by adopting both ladder [iisiiTiiiniiii] and 
cross-ladder approximations of the BS kernel [S]. No¬ 
tably, the extension to a bound fermionic system have 
been also undertaken [5]. It has to be recalled that 
numerical investigations of the homogeneous BSE has 
been performed also by considering the standard 4- 
dimensional variables m- 

The successful achievements for the homogeneous 
BSE encourage the extension of the Nakanishi inte¬ 
gral representation to the study of the inhomogeneous 
BSE, i.e. the integral equation that determines the scat¬ 
tering states. Our aim is to present a new applica¬ 
tion of our general approach [U, based on the so-called 
LF projection of the BS amplitude, i.e. the exact in¬ 


tegration on the minus component of the relative four- 
momentum that appears in the BS amplitude. After ap¬ 
plying this formally exact step to BSE, we have numer¬ 
ically investigated the zero-energy limit of the inhomo¬ 
geneous BSE, for a massive two-scalar system interact¬ 
ing through the exchange of a massive scalar, in ladder 
approximation. The calculated scattering lengths have 
been compared in great detail with the analogous ob¬ 
servables recently obtained Haul within a completely 
different framework. We have also compared our re¬ 
sults with the non relativistic scattering lengths, with 
the intent to yield a possible guidance for lowering the 
model dependence in the treatment of interacting fi¬ 
nal states, pertaining to relevant hadronic decay modes. 
Indeed, improving and widening our study could con¬ 
tribute to achieve an actual evaluation of the covari¬ 
ant off-shell T-matrix, that represents a key ingredient 
for describing, e.g., the heavy meson decay amplitudes, 
like in D —)• ATtttt processes p^fTH] . and could also have 
an impact in the development of final-state-interaction 
models, needed in the analysis of the CP violation in 
charmless three-body B decays HU. Moreover, we have 
properly analyzed the distorted part of the zero-energy 
wave function, putting in evidence the relation between 
a non smooth behavior of the Nakanishi weight function 
and the expected singularities of the LF 3D wave func¬ 
tion, like the one that brings the information relative 
to the global propagation of the interacting two-scalar 
system. Finally, the integral equation for the Nakan¬ 
ishi weight function obtained by applying the so-called 
uniqueness theorem m, is carefully analyzed for the 
general case of positive energy. Such an in-depth anal¬ 
ysis allows us to illustrate a surprising change in the 
width (from (— 00 , 00 ) to [0,oo)) of the support of the 
Nakanishi weight function with respect to its non com¬ 
pact variable, when the zero-energy limit is considered. 
Clarifying this feature allows us to put the forthcoming 
calculation of the phase-shifts on a sound basis, since 
their calculation requests a careful analysis from both 
the theoretical and numerical points of view, as it will 
be illustrated elsewhere HU- By concluding this Intro¬ 
duction, it could be useful to remind that developing 
genuine non perturbative descriptions of the scattering 
processes within the Minkowski space, possibly apply¬ 
ing formally exact frameworks, is an appealing goal, in 
view of attempts of extracting tiny, but fundamental 
signals once very accurate experimental data will be¬ 
come available. 

The paper is organized as follows. In Sec. we 
shortly introduce both the definitions and the general 
formalism, and we thoroughly discuss the problem of 
the support of the Nakanishi weight function, given its 
relevance for the zero-energy limit. Sec. illustrates 
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how to evaluate the scattering length from the Nakan- 
ishi weight function, in ladder approximation. In Sec. 
1 ^ the numerical studies of the scattering length are 
presented and compared with the existing calculations 
found in literature; moreover the scattering 3D LF wave 
function (indeed the distorted part) is analyzed. Finally, 
in Sec. the conclusions are drawn. 


defines Zi = —2kf /M and gets Zi = 2k~/M, since the 
incoming particles are on their-own mass shell: {p/2 ± 
ki)'^ = w?. Moreover, one has 1 > \zi\, since the in¬ 
coming particles have positive longitudinal momenta, 
i.e. p+/2±k+ > 0. In Eq. 0 . the following notations 
have been used: (i) cosp = kj^ • kij_, (ii) 7 = |kj^p and 
7 i = |kij_p, and (iii) = m? — M^/4. For the initial 
state one has 


2 The Nakanishi Integral Equations for 
scattering states 

In this Section, (i) we quickly recall the general formal¬ 
ism of Ref. [9] , for obtaining two integral equations that 
allows one to determine the Nakanishi weight function 
needed for scattering processes, and (ii) we demonstrate 
a relevant feature of the weight-function support, that 
it turns out to be very important also for numerically 
solving the inhomogeneous BSE. 

In our investigation, we considered an interacting 
system composed by two massive scalars that exchange 
a massive scalar. This is a generalization of the hon¬ 
orable Wick-Cutkosky model [20ll^ in two respects: (i) 
the interaction takes place through a massive-scalar ex¬ 
change and (ii) the scattering states is our focus. 


2.1 General formalism 


For scattering states, the incoming particles are on their- 
own mass-shell and we indicate their total and relative 
four-momenta with p and ki, respectively. By assum¬ 
ing that the inhomogeneous BS amplitude ‘P'^{k,p,ki) 
be expressed in terms of the Nakanishi weight-function 
, z', z” ■,K^,Zi), then one can write (cf Ref. [9]) 

<P^+Hk,p,h) = 

/ I ^1 poo 

dz' / dz" / d'j' 

-1 J-i J-00 


X " — 

[ 7 ' -I- m? 


g(+'>{i,z',z";K\z,) 

— ki^ — P ■ k z" — 2k ■ ki z' — ie] ^ 


= (27r)'^i5^^^(fc — ki) 



g(+)(7^z^z";K^^,) 

m -13 ’ 

[Vo - ie\ 

where the total four-momentum is p = {M, 0} and 
'Do = ^' + J + k'^ - k~{k~^ + - ^Ziz') 

-k'^^{z" + Ziz') ^-2z'cospy/ZfZfi ■ 

The power of the denominator is the same one adopted 
for describing a bound state (cf Refs. pilHI^ ITUlfTT] ). Ex¬ 
ploiting a standard formalism introduced in Ref. [3] , one 


(p/ 2 ±fci)^ = 


— + k+k, - 7 i = ( 1 - 2 ;J —- 7 ,, 

( 2 ) 


with necessarily (zi)^ < 1. Hence one gets 


= 4 


{w? -f 7 j) 

( 1 -^D 



< 0 , 


(3) 


To complete the generalities, we also give the expres¬ 
sion for the inhomogeneous BSE, without self-energy 
insertions and vertex corrections, in the present stage 
of our approach,. Then, one can write 


ki) = {2TT)^5^‘^\k - ki) 

+G^o^\k,p) y /C(fc,fc',p)^(+^(fc',p, fci), (4) 

where i /C is the interaction kernel (where the vertex 
corrections should appear), and is the free two- 

particle Green’s function given by 

(|+fc)2-m2 + ie (I-fc)2-TO2 + ie ■ ^ ’ 


It is worth noting that the bosonic symmetry of the 
BS amplitude, Eq. Q, when 1 —>■ 2, (i.e. p —>■ p, fc —>■ 
(—fc) and ki —)■ {—h)) has to be fulfilled, as in the case 
of bound states [10]. Therefore, the Nakanishi weight 
function must have the following property 


5 ^"^^ ( 7 ', z', z"; K^,Zi) = 5 ^+^( 7 ', z/-z"\ti^,-Zi). ( 6 ) 


Moreover, as shown in details in in [Appendix A[ one 
has 


( 7 ', = ±1, z") = ( 7 ', z', z" = ±1) = 0 . (7) 

As well-known (see e.g. Refs. E11I231SI), by projecting 
the BS amplitude onto the null-plane, i.e. integrating on 
k~, one exactly gets the 3D LF scattering wave func¬ 
tion that is proportional to the valence component 
^n- 2 /p cippearing in the Fock expansion of a two-scalar 

state, namely = \/2 V’i~!L 2 /p (given the normaliza¬ 
tions assumed in Refs. pini). The 3D LF scattering 
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wave function reads 

{z,'y,cosip;K^,Zi) = 


= p+- 




X {2Tr)^{k - ki) + i^dist {z,^,cosif,K^,Zi) (8) 

where k = {fc+,kx} and ipdist {z, J, cosp-, , Zi) is the 
distorted part of the 3D LF scattering wave function, 
that in the CM frame, where p'^ = P~ = M/2 and 
Pj_ = 0, reads 


4’dist (z, 7, cosp; K^,Zi) = 


(1-z^) 


dz' 


d'y 


,5(+)(y^zW^ 

[Vi - ief 


(9) 


with 

= y + 7 + + (1 - z^)k^ 

,,M 2 ^ ^ 

+Z Zj + 2cosp^77i) . 

In what follows, without loss of generality, we choose 
a head-on scattering process, namely a 2 -axis along the 
incoming three-momenta. In this case the variable 7 ^ 
is zero and therefore the dependence upon cosp disap¬ 
pears. As a matter of fact, the distorted wave function 
becomes 

fpdist {z,7]K^,Zi) = -- J dz' 


d'y' 


y+)(y,z',2;y,2,) 


[V 2 - ief 


with 


= 7 ' + 7 + z’^m^ + (1 - z'^)k^ + z' —z Zi 


and 


T 


( 10 ) 

( 11 ) 


Remarkably, tpdist (z, 7 ; Zi) displays a cut, originated 
by the free propagation of the two constituents, just as 
in the non relativistic case. In particular, the distorted 
part of the scattering wave function can be rearranged 
in order to make explicit the free propagation, obtain¬ 
ing (see details in [Appendix B ) 

I ( 2 \ ■ (1 - - 2 ^) 

Wdist{Z,7]K ,Zi)= I - - - 

X r dC dc' 

[k"^( 1 — z^j-I-m^z^-I-7 — lej y_i 7_i 

/ OO 

d7"^+(7",C'^C';^^^^) yi-in-IC'D 

-00 

(1 + z) d(C - z - C"z,) , 


-k 


(l + C'-C"^*) I?3(^,C',C")-*e 
(1 - z) d(z + C'z^ - C') 


(l-C' + C"^*) I?3(-^,-C',-C")-*eJ ’ 


( 12 ) 


where 

i?3(^,c',c") = «yi-n+mv+7 

(l-kz) 




(13) 


and ^+( 7 ', C") C'; -2^i) is the Nakanishi weight func¬ 

tion for the half-off-shell T-matrix (see Ref. i)- In par¬ 
ticular, the relation between the two Nakanishi weight 
functions is given by 


n(7' 


, z , z; K , Zj 


.) = 


da 


= ^| ^ I dc'g+{^,-A^^ 


f-l 


a a a 


x0{a-\z'\-\C'mi-a-\C-z-z'z.\). 


(14) 


with all the constrains on the variables explicitly writ¬ 
ten. Indeed, notice that the dependence upon z in the 
weight function 5 |^^.j( 7 ', z', z; k^, Zi) should be read as 
z-l-z'zi (cf Eq. ( 66 ) in [9] and Appendix B of the present 
paper). From Eq. the analogy with the non rela¬ 
tivistic case appears evident, once the familiar form of 
the global propagation is recognized. As a matter of 
fact, one has (1 > z^) 

1 


7 -|- z'^Tm? -|- (I — z'^)k‘^ — ie 

( 1 -zy 1 




(15) 


where Mq is the free mass of the two-body system given 

by 

-r t) 

^0 - ^ J’ ■ (16) 


^2_ 4 


It should be pointed out that the cut in ipdist is mir¬ 
rored in the integral equation determining the Nakan¬ 
ishi weight function, in particular in the part governed 


by the dynamics (see Eq. (18), below). It is useful to 


anticipate that the cut is canceled by the proper factor 
in the evaluation of the scattering amplitude. 

An issue of fundamental relevance related to ipdist in 


Eq. (10) (or to ipdist in Eq- (12)) is to determine the sup¬ 
port of the Nakanishi weight function ( 7 ', z', z; k^, Zi) 

(or equivalently 0 ^( 7 ', z', (/'; z^)) with respect to the 

non compact variable 7'. While the variable 7 = k']_ 
in Ipdist {z, 7', , Zi) is such that 7 S [ 0 , 00 ) and the 

same holds for 7 ' in the Nakanishi weight function when 
the bound state is discussed (see [ID]), in the case of 
a scattering state one has a different interval, namely 
7' G (— 00 , 00 ). Then, a question rises about the width 
of the support when —>- 0 “, i.e. the zero-energy limit 
which we are interested in. One should expect that the 
relevant support of 7' had to shrink in order to match 
the one pertaining to a bound state. 
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This can be accomplished if 


and F is 


^771^2^0- z; Zi) = 0 

for 7 ' < 0. Notably, this is what happens, as shown in 
detail in the following subsection. It should be pointed 
out that such a result is relevant for what follows, since 
we are going to consider the limit of a scattering state 
for 0 “, and one could be puzzled by the abrupt 

transition of the lower extremum for 7 ' from an un¬ 
bound value, for < —e, to a bound one, for the zero- 
energy limit. 


^ a(C'-^-^.C) , ( 1 -^)^ 

1, z] 7 ', C, C'; Zi) - ief (1 - C' + 

^ _ ojz + ZiC - CO _ 

p 4 iy, 1 : -z\ 7 ', C, -C'; -Zi) - ’ 

( 20 ) 

where 


2.2 The support of the Nakanishi weight function for 
the inhomogeneous BSE 


In order to address the support issue above introduced, 
let us consider the first meaningful approximation to 
Eq. Q, namely the approximation where the kernel 
i JC is substituted by its ladder contribution, given by 




2 

{k — — p? +ie' 


( 17 ) 


2^4(2/, 7, i, C, C'; Zi)=-i + -k K^(l - z^) 
+r{y,z,z^x,C n') + z{z,c,C;zi)—zzi (21) 

with 

r{y, z, Zi, c, c', 7') = 

X C) 7^ ^^) + ~ + + 7^1’ ) ( 22 ) 


First, one inserts the Nakanishi Ansatz for the BS am¬ 
plitude, Eq. 0 , in the ladder BSE. Then, one can 
perform the integration over k~ without any approxi¬ 
mation, and obtain the ladder inhomogeneous BSE pro¬ 
jected onto the null-plane, i.e. an integral equation that 
relates ipdist given by Eq. (10), to the dynamics dictated 
by the ladder kernel (see details in Ref. 0 ) . Namely, one 
gets 


,, , 9\ll){i,z',z-i^,zi) 




dz 


p2 - ie]^ 


9 




[7 -k z'^irF + (1 — z‘^)k^ — ie] - 
1 




/ OO ^1 

H J J ff (M) (7', C, C'; ) 

( 18 ) 


2 ( 47 r 

POO 

/ dy F{y,-f,z;j',C,C) 

Jo 

where is 

>V(^‘^^( 7 ,z;K^z 0 = 


1 


{z - Zi) 

9{z - Zi) 


^(l + z)-i^(l + z.) + ^-*. 


6 »(zj - z) 


^( 1 -Z)-^( 1 -Z.)+^-^ 


( 19 ) 


z{z,c,C;z,) 


(1 + ^) 

(1 + C ~ Zi() 


c, 


(23) 


and 


1\/t2 

AiC, c',7', ^^) «"(1 + C^) + 7 • (24) 


Because of the presence of the theta functions in Eq. 
( 20 ), one has 


1 >|C|>|^(±Z,C,±C',^*)I- 


(25) 


It should be pointed out that Eq. ( |l 8 | ) is relevant for the 
calculation of the phase shifts and, in the zero-energy 
limit, of the scattering lengths (cf Sec|^. 

After combining the global propagation with the de¬ 
nominator in ( 7 , z; k^, Zi) and repeating the same 

step for F (see Ref. [9]) one can apply the Nakanishi 
theorem on the uniqueness of the weight-function for 
an n-leg transition amplitude [12]. It should be recalled 
that the uniqueness theorem has been proven within a 
perturbative framework, while in the present context, 
a non perturbative one, the uniqueness is conjectured 
and numerically checked. Eventually, one gets a new in¬ 
tegral equation for the Nakanishi weight function, that 
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allows us to discuss the support issue, viz [9] 

9{-z') 5{'^ -'^a{z')) x 
^9{z - Zi) 9[l - z + z'{l - Zi)] 

■\-9{zi — z) 9 \\ -\- z -\- z'Zi)] j- + 

/ OO pi 

H j ^dc, J ^dC 


in the present work, cf Sec. leads to understand in 
detail the sharp transition of the support in 7 . 

For < 0 the support is (— 00 , 00 ), and one can 


split the integral equation (26) in two coupled integral 


equations: one is inhomogeneous, while the other is ho¬ 
mogeneous. To show this, let us introduce the following 
decomposition of the weight function 5 ( 21 ) (T'’ 


-k 


2 ( 47 r 

(1 + C - ZiQ 


+ ^'(- 7 ) ■ (33) 

Inserting such a decomposition in Eq. (p^ one gets 


(+) 


(1 ~ C' + ZiC) 


d'h, z', -z, -Zi; 7 ', C, -C', M ) 


where 

'faiz') = z'{2^^ - > 0, 

and h'{j”,z',z,Zi;j'X,C,9'^) is given by 

(1 + ^) 


(26) 

(27) 


h'(7",z',^,z,;7',C,C',5^) = 


(i + C'-^.C) 


'l9A 


pOO pi 

dy 6[z'- ^Z{zX,C'^z:i)] m ni ni 

Jo Jo + dy dC dC 

J^{X,yy;y',z,C,C,y;z,,y,^iy]} . (28) 


4l(id)(7>2:^z;K^z,) 9{-z') (I(7-7a(z')) X 

| 6 »(z - Zi) 0 [1 - z -I- z'(l - Zi)] 

+^(zi — z) 5 [1 -|- z -|- z^(l -|- Zi)] J- -|- 

X H'y,z',z,Zi;yX,C,yy 5p22d)(7'>C, C'; z:,) 


X 6 
with 

y, y,y\z, c, c, 7 '; z*, t^) = 

n ^ 

= 1 -i 


X iJ'(7,z',z,Zi;7',C,C',M^) 5lt(L)(7'>C.C';K^,z:,) 


(34) 


(1 + z) 


and 


(1 + C ~ ZiC) 

5M(C,C',7',^^)+5(5^+7') + M^ 


-^A (29) 


5it(L)(7,z:',z;K2,z,) = 9{-y 


2(47r)^ 


Notice that the inhomogeneous term vanishes both at 
z = ±1 and at z' = —1, as expected (cf Eq. 0 ). Indeed, 
for z = I, one has 

9{-z') | 6 »(I - Zi) 9 [z'(I - Zj)] 

+^(zi — 1) 9\y + z' {1 + Zi)] I, (30) 

that vanishes. For z = — 1, one gets 

9{-z') |6»(-1 - Zi) 9 [z'(l - Zi)] 

+9(zi -|- I) 0 [z^(l -l- Zi)]|, 

and again the theta functions produce a vanishing out¬ 
come. Finally, if z' = —1, the inhomogeneous term is 
vanishing, since 


dy / dc / dc 

'-Jo J-i J-i 

X H'y,z',z,zi;yx,c',yy 5p2L)(7'>C, C'; Zi) 


/ O pi 

dy / dc / dc 

-00 J-i J-i 


X H'y,z',z,zi;yx,c,yy 5it(L)( 7 ':C,C'; z*) 


(35) 


with 


(31) H'y, z', z, z,; y, C, c', = 




9{z - Zi) 0 [1 - z - (1 - Zi)] 

-|-0(zi - z) 0 [1 -I- z - (1 -I- Zi)] = 0. 


(32) 


(1 + C' - z.C) 

(1 - z)9{z -C + Zip 
(l-C' + z,C) 

xh'(7, z', -z, -zf, y, C, -C', d^) 


(36) 


The integral equation based on the uniqueness the¬ 
orem (that has been numerically verified for the bound 
states case in Ref. uni, and for the zero-energy limit 


If —)■ 0“, the off-shell kernel in the homogeneous in¬ 
tegral equation, namely the one with 7 < 0 and 7 ' > 0 , 
becomes vanishing and this leads to a system of uncou- 
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pled equations. As a matter of fact, one has for —)• 0 

6»(-7) B{i) 2;', z, Zi = 0; 7', C, C', = 

X ^ [-A(A, y, 7", z, C, C', Y; Y = 0, }, 


A=0 


= 0 , 


(37) 


since the delta function is always vanishing, given 7 < 0 
and 

(1 + z) ^j;M(C,C',7',Y=0) + y(/r2+y)+/,2' 

^ (1 + CO 

+C A = 


y 


= C 


(i + Y 
(i + CO 


+ 7'^ + ^2 


+C A > 0 . 
Then, for 


^n‘,{Ld) 




0 , Eq. becomes 
5^ 


= Zi = 0) = -; 


0(-7) 


2(47r)2 

[ di ( dC [ dC H'{'y,z,z,Zi;y,C:C,d'^) 

J-oo J-1 J-1 

^ (Y, C, C'; = Zi = 0). (38) 

The above homogeneous integral equation, valid in the 
zero-energy limit, is expected to have as a solution only 
(7^'*’^^^^)(7, z', z;/v2 = Zi = 0) = 0, given the freedom 
in choosing for scattering states. Let us recall that 
for the bound state case, where > 0 and 7 > 0, one 
gets a homogeneous integral equation and deals with an 
eigenvalue problem. In particular, one finds a discrete 
spectrum for g'^, once a value is assigned to and 
fi (see, e.g.. Ref. m. where the bound state case is 
discussed, within the present approach). 

It is also instructive to trace the behavior of the 
previous coupling term when approaches 0“. If 
is different from zero, than the delta function in Eq. 


(37) can give a finite contribution, since its argument 


can vanish. To achieve such a possibility, one must have 
(remind that 7' > 0) 

A/f^ 

a(c,c',Y,Y) = — c'^' + ^Yi + cY + Y < 0,(39) 

since the other terms, -|- 7' and A, always yield a 
positive contribution (A approaches zero from positive 
values). The above constraint leads to a volume of the 
integration in the space {7', C} (it is a hyperboloid), 
that shrinks to zero for > 0“, viz 




Y +YY + Y < -Y 


In conclusion, for scattering states in the limit —>■ 
0“, the corresponding Nakanishi weight function re¬ 
duces to the component z', z; = z^ = 0 ) 

and fulfills the following inhomogeneous integral equa¬ 
tion 


Y = Zi = 0) = Y di-z') 6{-f --fa{z')) 
X {6»(z)6» [1 - z -k z'] -k e{-z)e [1 -k z -k z']} 


2(47r) 


nOO pi 

, 0 ( 7 ) J dj' J dC 




X ffSL)(YY,C';Y = z, = o) (40) 

This is sharply different from the general case < 0 
given by Eqs. (34) and (35). 


3 The Scattering length 

In the CM frame, the differential cross section for the 
elastic scattering of two scalars can be written as follows 

ESI 

with Tf^ the invariant matrix element of the T-matrix, 
that is dimensionless (recall that in a (f>^ theory the 
coupling constant g has the dimension of a mass), and 
f‘^\s,0) the elastic scattering amplitude. It turns out 
that 

r'YY)= ( 42 ) 

where s = and cos 0 = kfki. To introduce the rela¬ 
tion with the phase shifts Si, let us expand the scatter¬ 
ing amplitude on the basis of the Legendre polynomials, 
Pi{cos 6 ), as follows 

f E(2^+ 1) ff ^Ycosd), (43) 

^ £ 

where the relative three-momentum is sj^ — w?^ 

or and the projected amplitudes are given 

by 

ff = sinSi . (44) 

Finally, in the zero-energy limit, only the amplitude 
with £ = 0 survives and one obtains 

fo (Jo ~ -a 4 , (45) 

where a is the s-wave scattering length. Therefore, in 
the zero-energy limit one gets 

lim /^YsY) = -a 

s—>- 4 m^ 


(46) 
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On the other hand, the scattering amplitude can be 
calculated through the BS amplitude as follows (see 
Ref. [3] for details) 


r'is,o) = - 


1 


^/s 8 tT k'^kf 


[7 


liin {k',p\GQ\p)\<P^-^'>;p,ki) = 


lim 


y/s 8 TT (7.z)-).(7/,z/) 

4 / / 2 

Vdist{z,^;K ,Zi 


+ z^rn? + (1 - z^)k^ - ie] 


(47) 


(l-z 2) 

where k' = {p[ — p'^j^, p'^ + p '2 = P (recall that p^ = 
= s). 

In ladder approximation and choosing 7^ = 0, from 
Eqs. (10) and (18) one gets 

/(l,)(s,0)= ^ [w(^‘^)(7/,z/;«^^.) 


+ 


2(47r) 


di' 


dC / dC 


dy F{y,-ff,Zf;-f”X,C) 5(Ld) ( 7 ": C, C'; , 

(48) 


where 


If z 


a = 


9 


m^lOTT 


±1 the free mass Mq 


00, and one can see 
that >V(“)(7,z = ±1 ;k^,Zi) —)■ 0, by taking into ac¬ 
count also the constraints generated by the theta func¬ 
tions. In general, the denominators in Eq. (19) do not 
vanish, since (i) only minus components of on-mass- 
shell particles are present there, and (ii) the momentum 
conservation law does not hold for those components 
(one can also explicitly check that the denominators 
do not have real roots). Moreover, since 7i = 0 then 
= m?IZi = {m? ^f)IZf. It is useful to introduce 
some kinematical relations relevant for describing the 
scattering process. In particular, the initial and final 
Cartesian three-momenta, kj and kj, has to be com¬ 
pleted giving the third components, viz 


kiz = ^ (fc+ - K) = -Zi Y 

= ^(*/-V) =-^f Y 


(49) 


Then, one can write down the relation between the scat¬ 
tering angle 0 and the LE variables z/ and Zi, given by 

k; • kf = ZiZfY~ = COS0, (50) 

where p-ki = p-kf = 0 has been used (those constraints 
are imposed by the on-mass-shellness of the particles in 
the elastic channel). It also follows that 


\k,\^ = \kf\^ = -n^ 


Finally, by exploiting the relation = —z^M^/4, that 
holds for 7i = 0, one gets 


Zf = Zi cosO , 


(51) 


For = {m? — s/4) —>• 0“, both Zi and z/ vanish 
(as well as 7/ = \kfi_\^), and one loses the depen¬ 
dence upon the scattering angle 0 in the scattering 
amplitude, namely one has a s-wave scattering, as it 
must be. The two functions, (7/, z/; z^) and 

C, C')> become 


lim W^^'^\‘^f,Zf = Zi cos9\k^,Z i) = 

d{zf - Zi) 


= lim 


{ 




e{zi - Zf) 




-} 


(Zf-Zi) 

1 


le 


(52) 


and 


lim F{y, 7/, z/ = z* cos6»; 7", C, C') = 


[j/2 _|_ y (-^2 _|_ ryO-j _|_ ^2 _ ' 


(53) 


Then, in the zero-energy limit Eq. (48) reduces to (see 
also Appendix C) 


lim /old(sY)= -a = 

s—>- 4 m^ 


= "■“ +2(4 


-| noo pi 


X dy 


0 " [ 2 /^A(C/ 7 ")+ 2 / (ti^+ 7 ")j 

(54) 


where the first term in the curly brackets leads to the 
scattering length in Born approximation, viz 


a 

oba = — m ^ 


(55) 


Moreover, 5 oLd(7"j CO is the Nakanishi weight function 
in the zero-energy limit. It can be obtained by solving 
two different integral equations as discussed in detail 
in [Appendix C[ where the whole matter is presented 
in a substantially simpler way than the one in Ref. [5] 
(notice that a mistyping in Eq. (103) of [S] has been 
fixed). In particular, the integral equation that links 
'ipdist to the dynamics governed by the BS kernel in 
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ladder approximation is given by 

Tj-- _ 

Jq [7 + 7" + - ze]2 

= £ /“ '>(-") 


M 


[7 + 7" + 


|0(^) 0 [1 - ^ - y7^2] ^ e[i + z- } 

poo p\ poo 

/ di / <ff^t^(7',C') / rf7" X 

Jo J -1 J-00 

- - 2 0 ( 7 ") X 


2(4t 


[7 + 7" + z^TV? — ie] 

d(C'-z) h',{i',z-iX',d^) 


(56) 


Notably, /ig(7", z; 7 ', 7) M^) is the proper kernel for a 
bound state with vanishing energy, as one can check in 
Ref. [TO] . 

The expression of /io(7"i 7^ Cj is given by (see 


details in Appendix C) 


h',ij",z-,j\C\d^) = 
= 0 


- l3o{z, C', 7 ', 7 ", - 2At\/+ 7 ' 

Bo{z,C,\l\i',d^) }_ 

-^o(C',7') 7 " 


(1 + z) 7y_ [y27lo(C', 7 ') + J/(m^ + 7 ') + 


(l + CO 


dy 


(1 + z) Jo " [y2^g(7^ Y) + y (^2 Y) ^ 2]2 ’ 

(57) 


with 


A(C',7')=C''™"+7' = C''™2+7' > 0 , 

So(aC',7',7",M^)=A^" + 7'-7"^^^ <0 , 

(1 + z) 

^ 0 ( 2 :, C , 7,7 ,M ) = ^o(2>C ,7,7 ) 

-VAo(c',y) > 0, 

_ 1 

2Ao(C',70 

X [-So(z,C',7',7",M^) ± ^o(2:,C',7',7",M^)] ■ (58) 


The zero-energy limit of Eq. (26), i.e. the integral 


function, reads (cf Ref. 

o2 


IS 


5^ti(7,^) = 72 ^(7)d[M^(l-kl)-7] 


2(47r) 


/ I poo 

^dC d7'<7^tk7',C') X 


(59) 


d(C'-x) /i77,z;7',C',/i^) 

+ ^^^ 7 ^ d(z-C') /io( 7 ,-z; 7 ',-C',Ai^) 

It should be pointed out that the presence of a non 
smooth behavior, like the discontinuity around 7 ^ 
/r^(l — |z|), is expected if one has to reproduce the sin¬ 
gular behavior of the distorted part of the scattering 
wave function (cf Eqs. (10) and fill). 

As illustrated in the next Sec.|^ we have taken profit 
of the general structure of the weight function suggested 


by Eq. (59) for obtaining numerical solutions of both 


Eq. (56) and Eq. (59), and eventually calculating the 


scattering lengths. 

It is worth noting that the scattering length given 

(+)/ 


by Eq. (54) represents a normalization for Q') , 

when fj, < 2m. As a matter of fact, from Eq. ( [5^ , one 
realizes that the inhomogeneous term is different from 
zero only for 

0<7<M^(l-kl) ■ 

Moreover, within the previous interval and y < 2m, the 
contribution to the kernel h'o that contains 


6 


Yj-Y^ - 7 ' - - 2Ai\/C'^w 2 -h 7 ' 


(l±z) 
disappears, since 

77;-^^ - 7' - - 2^y+ 7 ' < 

(l±z) v 

(1 ± zj 

< (1 ± C') - - 2Atm|C'| < y\C\ (.±d - 2m) j(60) 

The final step in the above expression is always negative 
when y < 2m. 

Therefore, for 0 < 7 < /i^(l — |z|) and /i < 2m, one 
has 

, , „2 2 rl roo 

= l^-V I d-,' 




/ dy - 

Jo [y2Mo(C',70+ + 7 ') 

= —167rm a , 


( 61 ) 


equation based on the uniqueness of the Nakanishi weight where Eq. (54) has been exploited in the last step. 
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4 Results 

In this Section, the numerical studies of both the scat¬ 
tering length and the distorted part of the 3D wave 
function are presented. First of all, let us illustrate our 
numerical method for solving the two integral equa¬ 


tions in (56) and (59). The main ingredient is the fol¬ 


lowing decomposition of the Nakanishi weight function 
that takes into account the singular behavior shown in 


Eq. (59), but also the result in Eq. (61), that holds for 


fi < 2m (this is always fulfilled for realistic cases) 


ffoLd(7, z)= ^ ei-t) + e{t) ^ ^ Atj Gi{z) Cj{t) , 

e=o j=o 

(62) 

where (i) t = 7 — /r^(l — |z|), (ii) the functions Ge{z) 
are given in terms of even Gegenbauer polynomials. 


N, Ng 


G 2 ^^'^\z) (recall that must be even in z) by 

G,(z) = 4 (1 - z^) E(5/2) 


d+) 


' (2£ + 5/2) (2£)! (5/2) 

'K^{ 2 ^ + ?>) 


21 


\z).. 


(63) 

and (iii) the functions Cj{t) are expressed in terms of 
the Laguerre polynomials, Lj{bt), by 

Cj{t) = Vb Lj{ht) . (64) 

The following orthonormality conditions are fulfilled 


J dz Ge{z) Gn{z) = 5i„ 

pOO 

/ dt Cj{t) Ct{t) = 

JQ 


= b 


dt e * Ljibt) Li{bt) = Sji. 


(65) 


In [Appendix D some details are given for illustrat¬ 


ing how Eq. (59) can be numerically solved by using 


the previous decomposition. In order to speed up the 
convergence, in the actual calculations the parameter 
b = 15.0/m^ has been adopted. The finite-range inte¬ 
grations (as those with respect to the variable z and the 
variable 7 when integrated up to /i^(l — |z|)) have been 
performed using a Gauss-Legendre quadrature rule. The 
infinite-range integrations, on the other hand, have been 
performed adopting a Gauss-Laguerre quadrature method. 


Finally, the convergence of the expansion given in Eq. (62) 


is very rapid, and adopting the values Nz = 10 and 
Ng = 24 well converged values have been obtained. All 
the results presented in this Section have been obtained 
for such a choice. Notice that at the end of the calcula¬ 
tion /3 resulted to be in agreement with the normaliza- 


In Tables Bi and 1^ the scattering lengths, Eq. 
(54), calculated by using the Nakanishi weight function 


obtained by solving both the integral equation (56), 


tion — 167ra shown in Eq. (61). 


apvsi and the integral equation (59), ajjpn, are shown 
for /i/m = 0.15, 0.5, 1.0 and values of the coupling 
constant a = /(IQTrm^), that range from a weak- 
interaction regime to a strong one. Moreover, for the 
sake of comparison, the results of Ref. m, acK, eval¬ 
uated within a totally different framework, based on a 
direct calculations of the half-off-shell scattering ampli¬ 
tude taking explicitly into account contributions from 
the singularities affecting the amplitude itself, are pre¬ 
sented in the second column. Eor reference, also the 
Born values of the scattering lengths are given in the 
fifth column. Erom the Tables, one can observe a very 
good agreement among all the three sets of numerical 
results, but some comments are in order: (i) the com¬ 
parison between apsv and ajjNi clearly confirms that 
the nniqueness of the Nakanishi weight function can be 
assumed with a very high degree of confidence, as we 
have quantitatively shown also for the bound-state case 
mm-, (ii) differences between acK [ 13 ] and our calcu¬ 
lations are present for /t = 0.15 m when the value of a 
approaches a value which corresponds to a bound state 
of zero-energy. In such a case, the scattering length di¬ 
verges (let us recall that, for the bound-state case, a 
is obtained as an eigenvalue of the homogeneous inte¬ 
gral equation, in ladder approximation), or there is a 
change of sign. Indeed, the above mentioned numerical 
differences do not represent a big issue (nonetheless it 
will numerically investigated elsewhere), given the com¬ 
pletely different theoretical frameworks adopted in Ref. 
m and in our work, and the well-known resonance be¬ 
havior of the scattering length, when a bound state is 
approaching a zero-energy scattering state. Einally, it 
is worth noting that the Born approximation gba rep¬ 
resents a quite good approximation only for small a 
(see also the following Eig. [^. Summarizing, the re¬ 
sults shown in Tables BH] and together with the 
calculations for the bound states unnKii], are a very 
strong evidence that the Nakanishi Ansatz, like the one 
for scattering states in Eq. ([^, represents a reliable 
tool for solving both homogeneous and inhomogeneous 
BSE’s in Minkowski space. 

In Eig. the scattering lengths for the above three 
values of g,/m are presented as a function of the ab¬ 
solute value of the scattering length in Born approx¬ 
imation {gbaI (see Eq. (55) and the last columns in 
Tables and §. Interestingly, in the same figure, it 
is also shown the comparison with the corresponding 
non relativistic scattering lengths, evaluated through a 
well-known expression (see e.g. mm), that exactly 
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Table 1 Comparison, for /r/m = 0.15, between the scatter¬ 
ing lengths (see Eq. | |54[ )) evaluated in Ref. m , acK (sec¬ 
ond column), and the ones, apsv (third column) and auNi 
(fourth column), calculated by adopting the Nakanishi weight 
function obtained from Eqs. ( |56[ l and ( |59[ ), respectively. All 
the calculations are in ladder approximation. The first col¬ 
umn contains the coupling constant a = /{IQ-wm?). Fi¬ 

nally, the fifth column shows the scattering length in Born 
approximation, Eq. ( |55| . The scattering lengths are in unit 
1/m.(*Private communication by J. Carbonell) 


a 

acK [IS] 

apsv 

auNi 

dBA 

0.01 

-0.460 

-0.459 

-0.459 

-0.444 

0.05 

-2.70 

-2.65 

-2.65 

-2.22 

0.10 

-6.92 

-6.66 

-6.66 

-4.44 

0.20 

-34.6 

-29.9 

-29.9 

-8.89 

0.30 

79.5 

105.9 

105.0 

-13.3 

0.40 

27.2 

28.0 

28.0 

-17.8 

0.50 

17.7* 

17.2 

17.2 

-22.2 

0.60 

12.8 

11.8 

11.8 

-26.7 

0.70 

8.66 

7.73 

7.72 

-31.1 

0.80 

3.73 

3.68 

3.68 

-35.6 

0.90 

-4.57 

-1.13 

-1.13 

-40.0 

1.00 

-28.1 

-7.89 

-7.89 

-44.4 

1.10 

900. 

-19.4 

-19.4 

-48.9 

1.50 

24.7 

66.9 

66.9 

-66.7 

2.00 

17.4 

23.2 

23.2 

-88.9 

2.50 

14.4 

12.3 

12.2 

-111.0 


Table 2 The same as in Tab. [ij but for the mass of the 
exchanged scalar g/m = 0.5. 


a 

a-CK US] 

O.FSV 

O.UNI 

a-BA 

0.01 

-0.0403 

-0.0403 

-0.0403 

-0.04 

0.05 

-0.209 

-0.209 

-0.209 

-0.20 

0.10 

-0.438 

-0.438 

-0.438 

-0.40 

0.20 

-0.971 

-0.971 

-0.971 

-0.80 

0.30 

-1.64 

-1.64 

-1.64 

-1.20 

0.40 

-2.50 

-2.50 

-2.50 

-1.60 

0.50 

-3.66 

-3.66 

-3.66 

-2.00 

0.60 

-5.34 

-5.34 

-5.34 

-2.60 

0.70 

-7.98 

-7.99 

-7.98 

-2.80 

0.80 

-12.8 

-12.8 

-12.8 

-3.20 

0.90 

-24.7 

-24.7 

-24.8 

-3.60 

1.00 

-103.0 

-103.2 

-103.0 

-4.00 

1.10 

62.0 

61.9 

61.8 

-4.40 

1.50 

11.0 

11.0 

11.0 

-6.00 

2.00 

6.34 

6.34 

6.34 

-8.00 

2.50 

4.54 

4.53 

4.53 

-10.00 


reproduce the second Born approximation, viz 
m asA 

m a = - - Tj—. -r 

1 + 


( 66 ) 


The chosen range of \aBA\ is [0,1.5], in unit of the in¬ 
verse mass m (the mass of the interacting scalars). Be¬ 
yond this interval, the scattering lengths can change the 
sign, as illustrated by the above Tables. Moreover, since 
m uba = —ct after fixing the value of /r/m one 

can follows the dependence of the scattering length on 


Table 3 The same as in Tab. [ij but for the mass of the 
exchanged scalar g/m = 1.0. 


OL 

acK US] 

apsv 

auNi 

dBA 

0.01 

-0.010 

-0.010 

-0.010 

-0.01 

0.05 

-0.051 

-0.051 

-0.051 

-0.05 

0.10 

-0.104 

-0.104 

-0.104 

-0.10 

0.20 

-0.217 

-0.217 

-0.217 

-0.20 

0.30 

-0.339 

-0.340 

-0.340 

-0.30 

0.40 

-0.474 

-0.474 

-0.474 

-0.40 

0.50 

-0.621 

-0.621 

-0.621 

-0.50 

0.60 

-0.784 

-0.784 

-0.784 

-0.60 

0.70 

-0.965 

-0.966 

-0.966 

-0.70 

0.80 

-1.17 

-1.17 

-1.17 

-0.80 

0.90 

-1.40 

-1.40 

-1.40 

-0.90 

1.00 

-1.66 

-1.66 

-1.66 

-1.00 

1.10 

-1.95 

-1.96 

-1.96 

-1.10 

1.50 

-3.79 

-3.80 

-3.81 

-1.50 

2.00 

-11.1 

-11.3 

-11.3 

-2.00 

2.50 

56.8 

51.7 

51.6 

-2.50 


the Yukawa coupling constant g^. In particular, from 
Fig. 0 one can see that for increasing values of and 
fj,/m the relativistic treatment in Minkowski space, be¬ 
comes more and more important, as expected, since the 
effect of the attractive interaction becomes more and 
more large. Notice that the scalar exchange in Eq. (17) 
leads to a non relativistic attractive Yukawa potential. 
Summarizing, modulo the adopted ladder approxima¬ 
tion, the comparison suggests that some care should be 
taken when one has to consider the effect of the interac¬ 
tion in the description of both hadronic scattering pro¬ 
cesses, even in the low-energy regime, and final states, 
that, e.g., are relevant for hadronic decays. 

In Fig. the Nakanishi weight functions for (i) 
g/m = 0.15, 0.5, 1.0, (ii) a = 0.1, 2.5, and (iii) z = 0, 
but running ^ jrr? ^ are shown. It should be pointed out 
that, for each value of g/m, the two values of the cou¬ 
pling constant a are representatives of a weak-interaction 
regime and a strong one. Moreover, since the Nakan¬ 
ishi weight functions obtained from Eq. (56) and Eq. 

(59) substantially coincide, only the calculations corre¬ 
sponding to Eq. (56) are shown. As mentioned at the 
beginning of this Section, the step-function behavior for 
small 7 has to be present, and the discontinuities are 
needed for obtaining the expected singularities in ipdist, 
like the one due to the global propagation. In Figj^ 
the transition from the weak regime to the strong one 
increases the discontinuous behavior, that for large fi 
become more and more smooth. Finally, recalling that 
for a bound state and g —>■ 0, i.e. the Wick-Cutkosky 
model mm, the Nakanishi weight function becomes 
proportional to i5(7), it is instructive to see the onset of 
such a behavior in the upper part of Figj^ 
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Fig. 1 The scattering lengths, calculated by using Eq. ( |54[ |, 
i.e. corresponding to solutions of the inhomogeneous BSE 
at zero-energy, vs |aByi| (Eq. ( |55[ )). Thick-solid line: a for 
/r/m = 1. Thick-dashed line: a for /r/m = 0.5. Thick-dotted 
line: a for /r/m = 0.15. The non relativistic scattering lengths, 
represented by the corresponding thin lines, have been cal¬ 
culated by using Eq. ( |66[ |. Notice that for fi/m = 0.15 the 
non relativistic calculation largely overlaps the relativistic one 
(thick-dotted line), and therefore it is indistinguishable. 


For 7 = 0 and \z\ ^ 1, only the first part of the 


decomposition in Eq. (62), i.e. 

|z|)), is dominant, and therefore trivial. 




In Fig. 1^ the same quantities as in Fig. but 
for 7/m^ =0.1 and running z, are also shown. As il¬ 
lustrated by the figure, the Nakanishi weight function 
acquaints a quite discontinues behaviour, as /i/m in¬ 
creases. 


Indeed, it is more profitable to present LF distri¬ 
butions, obtained from the distorted part of the zero- 
energy 3D scattering wave function. In analogy with 
the bound-state case (see Refs. [iniE]) , one can con¬ 
struct transverse and longitudinal LF momentum dis¬ 
tributions. In particular, one gets the following expres¬ 
sion for (z, 7; = z^ = 0) 


V'itf = Zi = 0) = 

4 Jo [7' + 7 + ■ 


(67) 


It should be noticed that inserting in Eq. (67) only the 
first part of the decomposition (621, one quickly reob¬ 
tains the singular behavior due to the global propaga¬ 


tion as shown in Eq. ( |12[ ), viz 
'IPdif 7; = Zi = 0) - /3 


(l-z2) 


1 


1 


7-b z^m^ /i2(l — |z|)-b 7 + z^m^ 

(i-z2)r 


= /3 


+ z'^w?) [g?(\ — \z\) + ^-\-z'^m?] 

( 68 ) 

Therefore, one has to expect singularities in the LF 
momentum distributions, that we would introduce in 
analogy with the ones for the bound states m- Let us 
emphasize, that only for the bound states they have a 
probabilistic interpretation. One could defines (i) the 
distorted transverse LF distribution 


V distil) = 


d£. 


2(27r)3 2 ^(1-0 


X dcj) (^- 7; = 0)]^ = ^ 

xrdz(i-z4rdy^ii^ 

7-1 7 o [7 + 7 + 


(69) 


and (ii) the longitudinal one, viz 
4 ^disti^) — 


(2^)3 2C(l-0 


j dk_L 7; = Zz = 0)]2 = 2 

Jo ^ Jo ^ [7' + 7 + ^W] 


(167r)2 


(70) 


with the fraction of longitudinal momentum given by 


^ ^ ^ (C+k+ 


(71) 


For the sake of presentation, it is useful to partially 
removing the singularities affecting the above distribu¬ 
tions. Therefore, in Fig. § 72 Vii) and |l-2^|3/2^(e) 
are shown. Figure illustrates the overall behavior of 
the LF distributions by varying the coupling a and the 
mass of the exchanged scalar /r, as in Figs. and 
It is worth noting the order-of-magnitude differences, 
when the coupling a is changed, but the typical fea¬ 
tures that one expects are still recognizable. A part the 
divergent behavior, already pointed out, that can be as¬ 
cribed to the global propagation, the transverse distri¬ 
bution shows the ultraviolet tail produced by the dom¬ 
inance of a single exchanged scalar, exactly as in the 
case of the corresponding distribution for bound states 
(see Ref. [H]). As to the longitudinal distributions, the 
expected peak at ^ = 1/2 or z = 0 is also seen. 
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|jym= 0.15 


)i/m= 0.15 




)i/m= 0.50 
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y/m^ 
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y Im 
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y/m^ 


Fig. 2 The Nakanishi weight function iz), in the zero-energy limit, vs 7 /m^, for /r/m = 0.15, 0.5, 1.0 and z = 0, 

Left panels: weak-interaction regime with a chosen value a = 0.1. Right panels: strong-interaction regime with a chosen value 
a = 2.5. 


The practical use of ijjdist is given by the evaluation 
of reactions that need a reliable treatment of the rela¬ 
tivistic effects, i.e. when the coupling constant becomes 
larger and larger. 


5 Conclusions 

In the present paper, our approach p iTOlITT] . based 
on the Nakanishi integral representation of the Bethe- 
Salpeter amplitude, is extended for the first time to 
the quantitative investigation of the zero-energy limit 
of the inhomogeneous Bethe-Salpeter Equation, in lad¬ 
der approximation, for an interacting system composed 
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|jym= 0.15 



z 


)i/m= 0.50 




Z 

Fig. 3 The same as in Fig. [^but for running 2 , and 7 = 0.1 

by two massive scalars that exchange a massive scalar. 
This achievement represents a non trivial task, that has 
allowed us to gain a sound confidence in the Nakan- 
ishi Ansatz, as an effective and workable tool for ob¬ 
taining actual solutions of the homogeneous and in¬ 
homogeneous BSE’s in Minkowski space. Indeed, the 
same approach that leads to a careful description of the 
bound states also yields a very accurate evaluation of 
the scattering length, as shown in Tables and by 


)r/m= 0.15 




z 


the quantitative comparisons with the same observable 
evaluated within a totally different framework, based 
on the direct calculation of the contributions from the 
singularities of the inhomogeneous BSE Haul]. 

As in the bound state case, we have performed the 
calculations by using the integral representation of the 
BS amplitude in terms of the Nakanishi weight func¬ 
tion, Eq. Q, that explicitly shows the analytic depen¬ 
dence of the BS amplitude upon the invariant kinemat- 
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|i/m= 0.15 


(i/m=0.15 








Fig. 4 The LF distributions obtained from the distorted part of the 3D LF scattering wave function (see Eq. ( |67[ l for 
/r/m = 0.15, 0.5, 1.0 and zero-energy limit. Left panels: transverse LF-distribution 7 ^ 'Pi'j) vs the 7 /m^ (cf the transverse 
LF-distribution expression in Eq.( |69[ )). Solid line: strong-interaction regime with a = 2.5; dotted line: weak-interaction reg ime 
with a = 0.1 Right panels: the same as in the left panel, but for |1 — 2^|®/^((>(^), (cf the longitudinal LF-distribution in Eq. 70 I 


ical scalars of the scattering process, under scrutiny. 
Then, by applying the LF projection onto the null-plane 
to the inhomogeneous BSE in Minkowski space, Eq. 
(without self-energy and vertex corrections), one is able 
to formally obtain the inhomogeneous integral equation 
for the Nakanishi weight function, that depends upon 
real variables. Its expression in ladder approximation 


is given by Eq. (181. Eventually, one can deduce an¬ 
other inhomogeneous integral equation for the Nakan¬ 
ishi weight function, Eq. (26), by assuming to be valid 
the uniqueness of the Nakanishi weight function, also 
in the non perturbative regime (recall that the theorem 
was demonstrated by Nakanishi m in a perturbative 
framework, but taking into account the whole set of 
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infinite diagrams contributing to a given n-leg ampli¬ 
tude). 

The numerical comparisons for the scattering lengths, 
obtained by using our Eqs ( |^ and ( [^ , and the cor¬ 
responding quantities calculated in Ref. |13j are shown 
in great detail in Tables □ i and 1^ It has to be em¬ 
phasized that the high accuracy reached by our calcu¬ 
lations is due to the new decomposition (62), suitable 
for obtaining the numerical solutions of the two inho¬ 
mogeneous integral equations, involving the Nakanishi 
weight function. The comparison with the non rela¬ 
tivistic scattering lengths (cf Fig. has illustrated the 
potential impact of a proper treatment of the relativis¬ 
tic effects in the investigation of hadronic scattering 
states, even in the low-energy regime. 

For the sake of completeness, the behavior of the 
Nakanishi weight functions, in the zero-energy limit, for 
weak- and strong interaction regimes have been shown 
in Figs. an di Those figures illustrate the non smooth 
behavior of the Nakanishi weight functions for certain 
ranges of the variables, that is an inheritance of the sin¬ 
gular behavior of the scattering states. Finally, we have 
defined LF momentum distributions, longitudinal and 
transverse ones, in analogy with the bound state case, 
(but without the probabilistic interpretation, entailed 
from the normalization of a bound state). Those distri¬ 
butions are shown in Fig. for the sake of illustration 
and reference purpose. It should be noticed that the 
transverse LF distributions show the expected ultravi¬ 
olet behavior, i.e. a power-like one, already found in the 
bound state case. 

In conclusion, the Nakanishi Ansatz for the BS am¬ 
plitude allows one to numerically solve in a very ac¬ 
curate way the inhomogeneous BSF, at least for the 
zero-energy limit. Such an outcome of our approach, 
together with the very nice results obtained for the 
bound-state case, strongly encourages to move to positive- 
energy scattering states, in order to evaluate the phase- 
shifts. If the phase-shifts evaluated within our approach 
(presented elsewhere [H]) will agree with the ones in 
literature [13], then the reliability of the Nakanishi Ansatz 
as a starting guess for obtaining exact solutions of BSE’s 
in Minkowski space could make a substantial step for¬ 
wards, confirming the great potentiality of this method, 
that can be applied to many other cases, changing di¬ 
mensions im, statistics, kernels, etc. 
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Appendix A: Boundary properties of the 
Nakanishi weight function ( 7 ', 2 ;', z") 


This Appendix is devoted to the analysis of the relation 
between (i) the Nakanishi weight function (7', z', z"\ 
that yields the integral representation of the distorted 
part of the 3D LF scattering wave function, and (ii) 
the weight function 0^(7",C)CO) fhat yields the inte¬ 
gral representation of the half-off-shell T-matrix (cf Eq. 
(58) in Ref. [9]). Notice that the dependence upon 
and Zi has been dropped for the sake of a light nota¬ 
tion. This analysis allows one to obtain the conditions 
fulfilled by ,z\z") when z' = ±1 and z" = ±1. 

It is worth noting that while for t/~*’(7", C) CO tti® 
constraint 0(1 — |C| — |C^I) holds, for the variables z and 
z' an analogous relation does not exist. 

The above mentioned relation between (7', z', z") 
and ^+(7", C, CO reads as follows (cf Eq. (63) in Ref. [2, 
where a factor of two is missing, as well as in Eq. (60), 
but it was not relevant for the formal discussion, since 
it can be reabsorbed in Q'^) 


g(+)( 7 ',z',z") = i2 f dai [ da 2 7 . - ^ ^ ^ 

Jo Jo (1 - ai - 02)'^ 
x 0 (l — ai — 02) 0(1 — oi — 02 — |z'| — |z" — oi -I- 02!) 

^C+f _ z" - \ 

\ (1 - Oi - 02) ’ (1 - Oi - 02) ’ (1 - Oi - 02) / 

(A.l) 


where 

z' = C (1 - Oi - 02) 

z" = C' (1 - 0i - 02) -I- Oi - 02 . (A.2) 

Notice that the constraints 0(1 — oi — 02) and 0(1 — 
ICI - IC'D = ^(1 - oi - 02 - |z'| - |z" - oi -k 02I) have 
been explicitly written, differently from Eq. ( 58 ) in Ref. 

m- 

In what follows it will be shown that the above theta 
functions lead to a vanishing Nakanishi weight function 
at |z'| = 1 or |z"| = 1. 

Given the presence of 0(1 —Oi — 02 — |z'| — |z" —Oi-I- 
02I) and 0 < Oi < 1, it is easily seen that for |z'| = 1 
one has 


5(+)(7',z' = ±l,z") = 0 


The same holds for |z"| = 1. Eirst of all, let us perform 
a change of variables, viz 


c 


Ol 


1 - (oi -k 02) , 

_ i-c + ^ 

2 


A — Ol — 02 

1 -C-A 


02 = 


(A.3) 


then Eq. (A.l) becomes 

X C iA + 


X Ctw, 


7 ' z' z"-A 


C’r c 


0(C-|z'|-|z"-A|). (A.4) 


For z" = 1, one gets z" — A > 0, since A G [—1,1]. 
Then 


0 (l-C-|A|) 0 (C-|z'|-l 

-bA) = 0(1 - C - 1^1) - (1 - C) - l^'l] = 0 

since A > (1—C)-l-|z'| > 0 and 1—C > |A|. For z" = —1, 
one has z" — A < 0, and then 

0(l-C-|A|)0[C-|z'| 

-(1 + A)] = 0(1 - C - 1^1) 0[-(i - C) - l^'l - ^] = 0 

since —A > (1 —C)-l-|z'| > 0 and 1 —C > |A|. Therefore, 
from the above results, one gets 

5(+)(7',^',^" = ±l) = 0 


Appendix B: The distorted part of the 3D LF 
scattering wave function 


In this Appendix, it will be shown how the expected 
global free propagation of the constituents can be fac¬ 
torized out in the expression of as in the non 

relativistic case. This result is relevant in two respects. 
On one side, it emphasizes the analogy with the non 
relativistic approach, and on the other side it allows 
one to understand the support of the Nakanishi weight 
function , z', z; 7^, z^), when 7' runs. 

In the CM frame (pj^ = 0 and = M), assuming 
without loss of generality a head-on scattering, i.e. 74 = 
0, the 3D LF scattering wave function, is given by [S] 

{z,-f]K^,Zi) =p+^^—^-^ J <P^+\k,p) = 

= ^ {2Tr)^- ki) -b V'dist(z,7; n^,Zi) , 

(B.5) 


where k = {fc+, kx} and ipdist is 


-ipdistiz,-/; K^,Zi) = 


(1-^^) r... 


dz 


d'y' 


'-1 


g(+)(7',z',z;«2,z,) 


[7' -b 7 -b z'^vn? -b (1 — z'^)k?' -b ^z z’zi — 


(B.6) 
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By using the Nakanishi weight function for the half-off- 
shell T-matrix, one gets the following expression [9] 

_l_ (1 — z^) f dk~ 


= p' 


27r 


(fc'^|Go(p)T(p)|fcf) = 


I (1 — z^) f dk 
= ' — 


(I — fc) — m'^ + ie J-i 


{^ + k) — m? + it 

pi pi poo 

dz’ / dC' / di 


/ — 1 J —oo 


^ z\ C'; Zj) 61(1 - \z'\ - Id) 

k^ + ^ — + ('p ■ k + z'2k ■ ki — + it 

Then, one can write 
ipdistiz,-/; K^.Zi) = 

“■"v: 

X e{i-\z'\-\c\) f ^ 


(B.7) 


= - p 


/ I pi poo 

dz' / dC / dj' 

-1 J-l J-oo 


1 


1 


(M/2 + k+) (M/2 - k+) (fc++ C'f + z'k+) 
1 

(f+fcF^TfTfc)^i^+^I7(W2+^ 

1 

"" (I - k)- - (I - k)on + ^e/(M /2 - k+) 

1 


I fe+(C'"+z'fc, )-'r-j'-K^+ie ’ 


(B. 8 ) 


where 

(^+k)- = 

4 M(l-z) 


ik++C'^+z'k+) 

2(to^ -I- 7 ) 

r 

2(to^ -I- 7 ) 


^2 M(l-hz) 


(B.9) 


with = —zM/2. Since kf = —k^ = —ZiM/2, one 
gets 

2 /* ^ poo 

tpdist{z,l-,K^,Zi) = dC dj' 

eii-\z'\-\c\) 


f dk~ 


iC-Z- z'Zi) 
1 


2^ f+fc--(H+fc)-„+i2e/[M/(l-z)] 
1 


f - A:- - (I - fc)o„ + z2e/[M(l + z)] 

1 

[k- - k^^ + i2t/[M{C - z - z'z^)] ’ 
with 

k- 2 

M(C' - z - z'Zi) 


(B.IO) 


M 2 


(C' + z'zi) +j + 'y' + , 


(B.ll) 


One has the following poles (recall that 1 > z > — 1) 


kh 


(| + k)or 



ku = 

km 



~ ^Na 


k)on + ^ + 

2 e _ 
M{C'^- z -z'Zi) ■ 


(B.12) 


In order to evaluate the analytic integration on k , one 
can consider the following two cases. 

If (' > z + z'zi, one can close the integration contour 
into the upper plane, taking the residue at ku, i.e. 


dk- 


1 


I 


1 


/_oo 27r [k - ku] [-k + ku] [k - kuu] 

1 

■ "* [M - (I - k)on - (i + k)on + It] 

1 


[f - (2 - k)on - kpf^ + it] 

- A 

(1-Z2) (I+Z) 


= -^^iC-z-z'z.) 


[k 2(1 — z^) + m'^z^ -|- 7 — ie] (I -|- C' — z'zi) 

I 


^(1_^2)+^2z2+7+ ^ ^ 


(B.13) 


If z -|- z'Zi > C,', one can close the integration contour 
into the lower plane, taking the residue at ku, i.e. 


dk- 


-00 27r [k - ku] [-k + ku] [k - kuu] 

- (C'-z-z'z.) 

, ( 1 -^) ( 1 -^^) 


(I — C'+ z'Zi) [k 2(I — z2)-I-m^z^-I -7 — ie] 
1 

(l-z)(M^zz'zi+-i''] 

^2(1 _ ^2) + ^2^2 + ^ ^ (Ll+z-z,) 


(B.14) 


where {(' — z — z'zi)t —)■ —e, since {(' — z — z'zi) < 0 . 
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Collecting all the above results, one gets the follow¬ 
ing expression for ipdist 

ipdistiz,-/; K^.Zi) = 

Al-z^) 1 


4 [k^(1 — z^) -\- m^z"^ -I- 7 — ze] 

I dz' I dC di ,z',Zi) x 
-1 J -1 J-oo 

(1 + ^) 

lil + C-z'z,) 

e{C-z-z'z.) ^(i-|z'|-|c'|) 

(l+z)( M^zz'zi+Y) 

^ 2(1 _ ^ 2 ) + ^ 2^2 + ^ ^ 
( 1 -^) 


— le 


+ 


(1 - C' + Z’z,) 

e{z + z,z'-C) 0(1-k'l-IC'D 


k '^{1 — z^) -I- m^z^ -I- 7 -I- 


(l-zjf^zz'zi+j') 

(B.15) 


One can reobtain the expression in Eq. (B. 6 ) by apply¬ 


ing the Feynman trick to Eq. (B.151. For instance, one 
has 

0(C' - Z - z.z') 


[k2(1 — z^) + m^z^ -|- 7 — ze] 

1 


9\ 9 9 (l + 2)(^Zz'Zi+7') 

K^(l - z^) -I- m^z^ - 1 - 7-1 - —- - «e 

= 0(C' — z — Ziz') f d^ X 

Jo 


(l+z)(Mj-zz'zi+'f') 

7 + m2z2 + k2(1 _ ^2) + ^ 

(l + C'-z'z.) 

a + z) 

_0(C' - Z - z,z') _ 

[7 -I- m'^z'^ + ^ 2(1 — z^) + a {^^zz'zi + 7') — ze]^ 

(B.16) 


f da 9 

\ (1 + ^) 1 
--^- ry 

Jo 

+ 

1 


with 


1 > 


(1 + ^) 


(1 + ^) 


(l-kC'-z'Zi) [1 +z+iC - Z-z'z^ 


since 0(C' — z — Ziz'). Inserting the above expression, 
together with the one containing 0 (z + Ziz' — C'), in 


distorted term 
'tpdist{z,r,Zi,K'^) = 

. (l-z2) /■! /■! r . „ /■' da 


/ I pL pOO pi 

dC dC di' / 

-1 J-1 J-oo JQ 




X 0(0-|C"|) 0(1-l-l-IC'D 

a 


[7 -I- ■mp-z'Z + k 2(1 _ 2:2) _|_ y/ _|_ ll|^2C"Zi — ze] 

{0 [(1 + z + C"z.) - 0(1 + C')] 0[a(C' - ^) - z.C]+ 
9 [{l-z- C"z,) - 0(1 - C'))] 0[-o(C' - ^) + ^.C"]} 

/‘a 


^ J—1 J-oo Jo ^ J—1 

x0(o-|z/D0(«-|C"D0(«-IC"|-M) 

Q+{T. C_ i/l 

^ ^ \ Q: ’ Q! 7 Q; / _ 

[7 -I- m?z'z + k 2(1 — z 2 ) _|_ y/ _|_ — ze] ^ 

{0 (1 -I- z -I- C"zi - a-y) 9{y - az - ZiC')+ 

0 (1 - z - C"zi -a + y] 9[-y + az + z^C)} , (B.17) 

where 7 " = 07 ' and C" = o z'. 

The theta functions between curly brackets single 
out the following integration regions 

— 1 — o -I- z -I- C'zi > y > az + ZiCJ' 

- az + ZiCJ' >y> -(1 - o) -I- z -I- C"zi 

The above intervals lead to the following constraint 
1 — a>y — z — C,”Zi > —(1 — o) 
namely 0(1 — a — \y — z — C''zi\) . Then one gets 

J ^ j j 


/-I 


a a a 


^ 0(a-|C"|-M)0(l-«-|z/-z-C"zD) 

[7 -I- wJz'z + «; 2(1 — z'Z) + 7" + ^z(”zi — ze] ^ 

(B.18) 

The above expression of ijjdist allows one to write the 
following relation between the Nakanishi weight func¬ 


tion ^+( 7 ', z', z;/t^, Zi), that appears in Eq. (B. 6 ), and 


namely the Nakanishi weight function involved in 
the description the half-off-shell T-matrix, 




5+( 7 ' 


.) = 


,Z ,Z]K , Zj 

da 7 ' z'y^ 

^ dy y^( — , —, -;k ,z,) 
/o J -1 a a a 


X 9(a - |z'| - |yD0(l - o - |y - z - z'zi|). (B.19) 


Eq. (B.151, one gets the following expression for the (A.4) by applying a suitable change of variables. 


Notice that Eq. (B.19) can be transformed into Eq. 
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Appendix C: Zero-energy limit 

The zero-energy limit of the relevant integral equations 
fulfilled by the Nakanishi weight function amounts to 
consider the case = 0, namely = 4 to^. This en¬ 
tails = Zi = 0 through = 4 (to^ -I- 7i)/(l ~ zf). 
In this Appendix, the integral equations obtained both 
without applying the uniqueness theorem [12j and by 
exploiting it, are obtained following a simpler proce¬ 
dure than the one adopted in Ref. [5] (notice that a 
mistyping present in Eq. (103) of 0 has been fixed in 
this Appendix, as explained in what follows). 

The Nakanishi integral equation, involving ipdisti for 

< 0 (see i) is given by 


/ I POO 

dz' J di' 5(l])(7", z', z\ "fi, Zi) 

1 

[7 -I- 7 " -I- z'^w? -b (1 — z'^)k^ + z'^zzi — ie] 

/ I POO 

dz' / dj" 0{-z')S{y'-Uz')) 

-1 J-OO 


[7 + 7 " -b (1 — z^)k^ -b -b — 
|6»(z - z^) 9[1- z + z'(l - Zi)] 

+ 6 [zi — z) 0 [1 -b Z -b Z^(l -b Zi)] J- 


z' ^zzi — ie\ 


9 

'2(47r)2 


/ OO p1 

dj” / dz' 
-00 J —1 


[7 -b 7" -b z^m^ -b K^(l — z^) -b z'^zzi — ie] 

/ OO p\ p\ 

H J d( J dc'5(Ld)(T^C,C';«:^^i) 

r (1 + ^) 

4(1 + C'-Z,C) 

X e{C - z - ZjC) /i'( 7 ",z',z,Zi; 7 ',C,C',Ai^) 

, (1-^) 

(1 ~ C + ZiC) 

X e{z-c + z^c) h'{y',z',-z,-zi]-f',c,-C:d'^) , 

(C. 20 ) 


with 


h'(7",z',z,z,;7',C,C',A^") = 


(1 + ^) 

(i + C'-^*C) 


d [J-(A, y, C, 7 ", C, C', l';z„^, y^)] . (C.21) 


where 

•^(A, y, 7", 2, C, C', 7'; z^, /r^) = 

c (1 + ^) 


= 7 - e 


(1 + C' ~ Zi() 


y^AiC,C,i\^^) + yid^ + i) + d 

y 


- CA (C.22) 


For = 0, it follows that 7 i = = 0 since 

^ \2 jvf2 


2 2 
n = rn — 




Then, taking into account that 3(7) (O'’ ~ 

has support only for positive 7 , one can write (cf Eq. 


( |C.20[ ) and subsec. \2.2\ 

J-1 Jo ^ [y + y' + ^2m2-ie]2 


= 9 


di' Cdz' ^(-')d(7" + .V) 


1-1 [ 7 -b 7 "-b z^m^ — ie]^ 


0 {z) 0{l-z + z') + 0{-z) 0{l + z + z') 


2 ( 47 r) 


di j ^dC 


X l\dC 3j7)(7',C,C';«" = 0) 


/-I 
poo 

X / dj‘ 


1 


dz' 


X 


[7 + 7 " -b z^m^ — ie] J-i 
i^0(C'-z) Z'(7 ",z',z;7',C',3^) 


where the kernel Z' is given by 

Z'(7",z',z;7',C'>7)= 

r*l 


(C.23) 


(1 + CO 

(Jx/ - Cdo(y,^,C'>70 - CA] 


xd 

with 


h^i+co^JjA=o 


(C.24) 


do(3,z,C'>70 = 


(1 + z) 1 


(C.25) 


(1 + CO y 

X Aio(C070 + y(M^ + 7 O + 7} 

and 

A/f2 

Alo(C',70=C''^+y= C'V+7' > 0, (C.26) 

Notice that 7 ' is positive and therefore also Fq has to 
be positive. Einally, 7 " in Eq. (C.24) has to be positive 
for getting a non vanishing Z', z', z; 7 ', /r^). 
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Performing (i) the integration on z' in both sides of with (see also Eq. (C.26)) 


Eq. (C.23) (recall that 1 > |^C(1 ±z)/(l±(^')|) and (ii) 


the integration on C in the rhs, one gets 




/g [7 + 7" + 

.2 roo 


9 


= ^ di 


n 


M J -00 [7 + 7 " + — ie]^ 

X {^(^) z-i' l\F\ 

+ 0 (_^) 0 [1 + ^ - J ^ 

^00 

H J ^dC gl^ldh'X') 


2(4t 


/ OO 

dq" 0 ( 7 ") 

-00 


1 


[7 + 7 " + z^TO^ — le]' 
1(11 ' /•! 2\ 


0(7 - 2 r) d^q", z; 7 ', C', 

where 

1 . 


(C.27) 


2 . 


50 ^( 7 ",^) = J ^dz' 9 [t\){l”,z',z;K^ = 0).(C.28) 

J dz' 6{—z') S{'y" + = 


3. 


= ^0(7")d(d^ - 7 ") • 

d 

u' i " ' /■' 2^ (1 + ^) d f 

hoil ,.;7,C,d)=(j^ 

[ d^d[7"-CEo(y,z,C',7')-CA] 
do 


(C.29) 


dy 


A=0 


(C.30) 

From Ref. uni, one recognizes that dg( 7 ", z;^'X\fj,‘^) is 
the suitable kernel for a bound state with zero energy. 
Therefore one can write 

do(7",x;7',C',d^) = 


= 9 


- Bo{z,C,i',i”,9^) - + 7 ' 


6o(z,C',7',7",A^^) 


_ ^ _i_ (1 + CO 

-dlo(C',7') A\o(2:,C',7',7",d^) 7 " (1 + 2 ) 


fy+ 

X / dy 


V- [yMo(C',70+ + 7 ')+ d^] ^ 


(1 + 0 


dy 


(1 + ^) Jo ^ [y'^AoiC,!')+ 9 ( 9 '^+ 1 ') + ' 

(C.31) 


Sg(z,C',7',7",0)=0+7'-7"^^^ <0 , 

(1 + z) 

Z\^(z,07',7",0) = 

= S2(z,07',7",0)-40dlo(C',7') > 0 , 

1 

2A(C',7') 

X [-eo(z,C',7',7",M')±^o(^,C',7',7",0)] (C.32) 

Notice that in the inhomogeneous term in Eq. ( |C.27 ) 
the factor 9{fA — 7 ") has been dropped, given the pres¬ 
ence of the step functions 0 [l ± z — q"//!^]. In Eq. 
(103) of Ref. ig the step function 9{fA — 7 ") has been 
accidentally overlooked. 

In conclusion, by applying the Nakanishi theorem 
on the uniqneness of the weight function [ 12 ) . one has 
the following integral equation 

0OLd(7, ^) = 77 ^( 7 ) d [0(1 - kl) - 7 ] 

„2 poo pi 

I di'J_^dA X 

0(z'-z) 0(7,z;y,z',0) 

+^ 1 +:^ 0(7, 7 ,- 2 :', 0 ) 


(C.33) 


Appendix D: An effective decomposition of 

9oid(7,^) 

In this Appendix, the decomposition of 59 ^^ ( 7 , z) shown 


in Eq. (62) is applied to the simple case of Eq. (59), 


based on the Nakanishi uniqueness theorem [12], in or¬ 
der to give the explicit representation of the numerical 
system to be solved. 


Inserting the decomposition (62) 


N. Ng 


0oLd(7,^) = P 9{-t)+9{t) ^^0 *^^0) 0(0, 

t =0 j =0 

(D.34) 

with t = 7 — — 71), in Eq. (59), given by (notice 


that in the following expression, the symmetry proper¬ 
ties of both the weight function and the kernel dg are 
exploited), 

0OLd(7, ^) = 77 ^( 7 ) d [0(1 - \z\) - 7 ] 

poo pi 

0 ( 7 ) / dq' / dz' 5^td(7', z') x 
Jo J -1 


{Any 

d{z' - z) do( 7 , z; 7 ', z', 0 ) 

(1 + +) 


(D.35) 
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one can quickly obtain the following coupled system 

p 9{-t) = ^ 0 ( 7 ) e{-t) - 0 ( 7 ) e{-t) p 




( 47 r) 




X / dz' I di 

J-1 ^0 

(1 + z') 


9 


N, N. 


(47r)^ 


%) EE ^ 




dz' 


e=o j=o 

d7' G'£(z') £j [7' - M^(l - |z'|)] X 


^*(2' - z) /io{ 7 , 7', z', ii^) , 

(1 + z') 
and 

EE = 

j=0 

.2 /.I /•M'li-b'l) 


(D.36) 


9 


(47r)^ 


/■i /-M (i-b I) 

0(7) d(t) P / dz' / 1^7' X 

J-l Jo 


^'(2' - z) h'oi-f, z; 7', z', + 

(1 + z'j 


{Any 


Ng 

^( 7 ) EE 

£^o i=o 




'-1 


/ d7' G£(z') Cj [7' - ^ 2(1 - |z'|)] X 

V(i-b'l) 

^(2' - z) /io( 7 ,2:; 7', z', yP). 

(1 + z'j 


For 7 = 0 and z = 0, Eq. (D.36) reduces to 


N, Ng 


P - —^ + P -f/3,/3 + E E Mj 

^ £=0 j=0 


(D.37) 


(D.38) 


where 


IfiR = — 


( 47 r) 


/ dz' 

'0 Jo 


l-t.pl-z') 


d'y' 


jyyrypj ^0(7 = o,z = 0;7',z',^ ) 


(D.39) 


and 


Ip.ej = - 


(47r)2 


pi pOO 

X dz' dj' Gi{z') £j [7' - pP{l - z'j] 

Jo J /i^ (1 —2') 


7—— /i'q(7 = 0, z = 0; 7', z', n^) 
(1 + z'j 


(D.40j 


with (cf Eq. ill)) 

K^l = 0 , z = 0 ; 7', z', /i^j = -(1 + z'j 


X dy 


0 (7'+ m^z'^ j + 1/(7'+ ^2) _|_ ^2j 


2 ■ 


(D.41) 


A matrix representation can be obtained for Eq. (D.37) 
by multiplying both sides by Gi> (zj Cj> {t) and integrat¬ 
ing, viz 

N, Ng 

Ai'j' = P +EE Aij (D.42j 


£=0 j=0 


where (cf Eq. 




/ I pQO pi 

dz / dt Gipz) Cjpt) / dz' 

-1 Jo Jz 

n-I-zj 

(D.43j 

0 0 + z') 


and 


Ii'3'd3= 


/ I poo 

dz / dt Gif{z) Cj'{t) 

-1 Jo 

pi poo 

/ dz' / dy' Ge(z') x 

Jz Jn^(l-\z'\) 

U3 b' - 9^0-\z'\)] K{l^z]-i',z',y?). 


(D.44j 
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